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The magnetic susceptibility measured in neutron scattering experiments in underdoped 
YBa2Cu307- y is interpreted based on the self-consistent solution of the t-J model of a Cu-O plane. 
The calculations reproduce correctly the frequency and momentum dependencies of the suscepti- 
bility and its variation with doping and temperature in the normal and superconducting states. 
This allows us to interpret the maximum in the frequency dependence - the resonance peak - as 
a manifestation of the excitation branch of localized Cu spins and to relate the frequency of the 
maximum to the size of the spin gap. The low-frequency shoulder well resolved in the susceptibility 
of superconducting crystals is connected with a pronounced maximum in the damping of the spin 
excitations. This maximum is caused by intense quasiparticle peaks in the hole spectral function 
for momenta near the Fermi surface and by the nesting. 



PACS numbers: 71.10.Fd, 74.25. Ha 



I. INTRODUCTION 

Inelastic neutron scattering experiments give im- 
portant information on the anomalous properties of 
high-T c superconductors. Among the results obtained 
with this experimental method is the detailed informa- 
tion on the magnetic susceptibility in YBa 2 Cu307_ y 
measured in wide ranges of hole concentrations and 
temperatures^*^ These measurements revealed the sharp 
magnetic collective mode called the resonance peak. The 
peak first observed^ in the superconducting state of 
YBa 2 Cu307 was later also detected in the underdoped 
compounds, both in the superconducting and normal 
states^*! Recently the resonance peak was also observed 
in Bi 2 Sr 2 CaCu 2 8+5 and Tl 2 Ba 2 Cu0 6+5 £ 

Theoretical works devoted to the resonance peak were 
mainly concentrated at the overdoped region where the 
peak is observed in the superconducting state and dis- 
appears in the normal stateib 2 i 3 i 4 i 5 i 6 In Refs. iTlSl lil an 
interpretation of the resonance peak based on the itiner- 
ant magnetism approach was proposed. This approach 
which uses the Lindhard function for the bare susceptibil- 
ity xo(kw) and the random phase approximation relates 
the appearance of the peak to the disappearance or con- 
siderable decrease of Imxo(Q^) hi the frequency range 
uj < 2A S with the opening of the d-wave^S or s-wave- 
superconducting gaps A s . Here Q = (it, n) is the antifer- 
romagnetic wave vector. In this frequency range the peak 
arises due to the logarithmic divergence in Rexo(Q^) 
which originates from the jumpi^ in Imxo(Qw) or due 
to the nesting of the bonding and antibonding Fermi 
surfaces in the two-layer crystal^ In the normal state, 
when the damping increases, the peak is smeared out. 
For the approaches of Refs. [8|0j the two-layer structure 
of YBa 2 Cu307_ a is of crucial importance for the ap- 
pearance of the resonance peak. However, the recent 
observation^ of the peak in single-layer Tl 2 Ba 2 CuOe+5 



indicates that an interaction between closely spaced Cu- 
O layers is not the necessary condition. A single-layer 
system described by a modified s-f Hamiltonian was con- 
sidered in Ref. In that work the resonance peak ap- 
pears also due to the vanishing damping of spin exci- 
tations and a real part of the spin excitation frequency 
stems from the fermion bubble. A qualitatively different 
approach was suggested in Ref. ll] where the existence 
of a well-defined branch of spin excitations which exists 
even in the absence of mobile carriers was postulated near 
the M point (k = Q) of the Brillouin zone. In this sce- 
nario the resonance peak is related to the excitation with 
k = Q of this branch. As in Ref. 0130 here the peak is 
visible in the superconducting state due to the absence 
of the damping for to < 2A S and is smeared out in the 
normal state. 

This latter scenario seems to correspond most ade- 
quately to available experimental data on the resonance 
peak. As mentioned, it is observed also in underdoped 
YBa 2 Cu3C>7_ a and in the superconducting state the peak 
varies continuously on passing from the underdoped to 
overdoped region. Moreover, in the underdoped region 
the peak is also observed in the normal state and its 
frequency is nearly the same as in the superconducting 
state. Therefore it is reasonable to search for a unified 
explanation for the peak which is applicable both for un- 
derdoped and overdoped regions, and - in the former 
region - for the normal and superconducting states. The 
existence of the excitation branch of localized Cu spins is 
well established for the underdoped regioniSiiii In this re- 
gion it is quite reasonable to connect the resonance peak 
with these excitations. 

In this paper we use the two-dimensional t-J model 
to which the realistic three-band Hubbard model of the 
Cu0 2 planes can be mapped in the case of a strong on- 
site Coulomb repulsioniH For the underdoped region the 
self-consistent solution of the t-J model was obtained in 
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Ref. [lj with the use of Mori's projection operator tech- 
nique. We employ this result for the calculation of the 
magnetic susceptibility in the normal and superconduct- 
ing states. The calculations reproduce correctly the fre- 
quency and momentum dependencies of the susceptibil- 
ity and its evolution with doping and temperature in 
YBa2Cu307_j / . This allows us to relate the resonance 
peak with the excitation branch of the localized Cu spins 
and to identify the frequency of the peak with the size 
of the spin gap at the M point. In the underdoped re- 
gion, with increasing doping the peak frequency grows 
with the gap size and the peak intensity decreases, in 
agreement with experimental observations. 1 '^ Moreover, 
the low-frequency shoulder observed^ in the suscepti- 
bility of superconducting crystals can be connected with 
a pronounced maximum which we find in the damping 
of the spin excitations. This maximum is caused by in- 
tense quasiparticle peaks in the hole spectral function for 
momenta near the Fermi surface and by the nesting. 



II. THE HOLE GREEN'S FUNCTION IN THE 
SUPERCONDUCTING STATE 

The Hamiltonian of the 2D t-J model readsi^ 



- V I (s z s z + s +1 



where a na = |ner)(nO| is the hole annihilation operator, n 
and m label sites of the square lattice, a — ±1 is the spin 
projection, |ner) and |n0) are site states corresponding to 
the absence and presence of a hole on the site. These 
states are linear combinations of the products of the 
3d x 2 _yi copper and 2p a oxygen orbitals of the extended 
Hubbard model. 13 In this work we take into account near- 
est neighbor interactions only, t nm = —t ^ a <5 n ,m+a and 
Jnm = J J2 a ^n,m+a where the four vectors a connect 
nearest neighbor sites. The spin-i operators can be writ- 
ten as sjj = \ J2 a 0"|ncr)(ncr| and s£ = |ncr)(n, —cr\. 

To investigate the magnetic susceptibility of this 
model the hole G(kt) = -i9(t)({A kr7 (t), A^}) and spin 
D(kt) = — i6*(i)([s k (i), sl k ]} Green's functions have to be 
calculated. Here supposing the singlet superconducting 
pairing we introduced the Nambu spinor, 
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thus G is a 2x2 matrix (here and below matrices and 
vectors are designated by boldface letters). In the above 
formulas the angular brackets denote averaging over the 
grand canonical ensemble, and 



a kff = iV- 1/2 V e- lkn 6 



ak<r(i) = exp(i7Y£)ak(T exp(— iHt), N is the number of 
sites, H = H — /i^ n I n , ii is the chemical potential, 
X n = |nO)(nQ|. 

To derive the self-energy equation for the matrix 
Green's function G the continued fraction representation 
for Green's function and the recursive equations for their 
elements from Ref. have to be generalized for the case 
of matrices. Such generalization reads 

Rn(w) = M-E„-R„ + i(w)F„]^ 1 , n = 0,1,2... (2) 

where I is a 2x2 unit matrix, the matrices E„ and F„ 
are calculated from the recursive equations 

[A„, H] = E„A„ + A, i+ i + F„_iA„_i, 
E^HUA^ffLAlJMA^At})- 1 , (3) 

F„ = ({A n+ i, A| l+1 })({A„, A^}} -1 . 

Here F_i = and Ao = Ako- for the case of the function 
G. As follows from Eq. Q), the two-component oper- 
ators A„ constructed in this recursive procedure form 
an orthogonal basis. For the anticommutator Green's 
function G the inner product of two arbitrary oper- 
ators A and B is defined as ({A, B}} and the or- 
thogonality means ({A„,A m }} = 5 nm ({A n , A„}). In 
Eq. 0, R„(w) = -if™dtexp{icut)Tl n {t), R„(i) = 
({A nt , At})({A n , Ajj}) -1 where the time dependencies 
are determined by the equation 

,— A nt = Y[(l P k )[A nt ,H], A n>t =o = A n 



k=Q 



4 = iV-^e^X, 



with the definition P„Q = ({Q, A+ }}({A n , A+ })~ 1 A fl 
of the projection operator P n that projects an arbitrary 
two-component operator Q on A n (for a more detailed 
discussion of these equations see Ref. 14). 

From the above definitions it follows for the Fourier 
transformation of G(ki) 

G(kw) = [wl - E - RiF ]- 1 <{A 0) Aj}), (4) 

where ({A , Ag}) = ipl, ip = i(l + x), x = (X n ) is the 
hole concentration, 

<TKlip-\ZJ^-U) -(e k -/0 )> 

(5) 

e k = -(4t<^+6iCi^- 1 +3JF 1 i^- 1 )7 k , // = f i+4tF 1 ip- 1 + 
SJCiip^ 1 , 7k = 4 5^ a exp(ika). The nearest-neighbor 
correlations C\ = (s+ 1 Sn| a ), Fi = (4 ff tt n +a,j), K\ = 
(7(a n ,a„+ a _j) and the hole concentration x can be ex- 
pressed in terms of the components of the hole and spin 
Green's functions: 

1 f°° 
x = — / dunp{io)A(\ujj), 

1 f°° 
F\ = T^X! 7k / dujn F (uj)A(kuj), 

I. J— oo 
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r <kjcoth(^)B(kLo), (6) 

i„ o 
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7k 



— 71^ (tL>)] 



x[L(kwcr) + iM(kw(7)], 



where 



L(kwCT) = -Im[Gi 2 (kcjCT) + G 2 i(kwo-)]/(27r), 
M(kwo-) = Re[Gi 2 (kwcr) - G 2 i(kwcr)]/(27r), 

A(kw) = -ImGii(kw)/7r, and S(kw) = -Im£)(kw)/7r 
are the hole and spin spectral functions, uf(uj) — 
[exp(a;/T) + and T is the temperature (the functions 
v4(ktj) and _B(kw) do not depend on a). In the deriva- 
tion of Eq. © for C\ we have taken into account that 
the approximation used retains the rotation symmetry of 
spin components 14 and therefore G\ = 2(s^s^ +a ). 

As follows from Eq. (JSJ), in the t- J model the supercon- 
ducting gap has an s-wave component if K\ ^ 0. How- 
ever, in the considered case this component is small in 
comparison with the d-wave component introduced be- 
low and will be neglected. 

From the definition of the hole self-energy X = RiFo 
we find that S 22 (kw) = EJ x (k, —lo) where it was taken 
into account that these components of X do not depend 
on a and are invariant under the inversion of k. For 
Sn(kw) the following expression obtained in Ref . 1 1 ll can 
be used: 



Im En(koj) 



167rf 2 x ^ 



duj' 



7k-k' + 7k 



7k) 



1+7V 



ReE u (kw) = V r — 

J-oo T 



sgn(w')(7k-k' ,, . 

V 1 ~ 7k' 

x [n B (-w') +n F {u-u')\ (7) 
x A(k-k',w-a/)5(k'a/), 
cko' Im£ii(ka/) 



where ris(cj) = [exp(w/T) — 1] and V indicates Cau- 
chy's principal value. 

Assuming the <i-wave superconducting pairing, for the 
anomalous self-energies we set 

Ei 2 (kwcr) = £i 2 (kwcr) = <rA"[cos(fc x ) - cos(fc a )]/2, (8) 

with the superconducting gap A s . For such anomalous 
self-energies M (kwrr) = 0. 



III. THE SPIN GREEN'S FUNCTION 

In Ref. ^| we have noticed that the approximation 
used there leads to an underestimation of the imaginary 
part of the magnetic susceptibility at low frequencies. 
To avoid this drawback in the present work we shall not 



split the spin self-energy into the hole and spin parts, 
but rather continue the calculation of the terms of the 
continued fraction using the entire Hamiltonian JJJ. 

The spin Green's function is calculated from the rela- 
tion 



D(M = ^((4|s-k))*-(4,sik), 



where 



Jo 

and Kubo's relaxation function 



(9) 



(10) 



((«)«= dte^ dt'([sUt'),st k ]) (11) 
Jo Jt 

can be represented by a continued fraction which is simi- 
lar to the scalar form of Eq. (J2J) . The elements E n and F n 
of this function are calculated from a recursive procedure 
which is similar to the scalar form of Eq. where, how- 
ever, mean values of anticommutators have to be sub- 
stituted by inner products of the type of Eq. (fTU|) (sec 
Ref. [H. 

From this definition we find for the starting operator 
Aq = s k of this recursive procedure 

E = (isl S z _ k )(sl,s\)- 1 =0, 



where 



"k>#]> 



, ... F _ 4(l- 7 k)(J|Ci|+*F 1 ) 

Using these elements of the continued fraction represen- 
tation of ((s k |sl k ))a>) Eq. © can be rewritten as 



D(kw) = A{1 ;*l J ™ + !?\ 



lu 2 - wll(kw) — UJ 2 . 



(12) 



where w k = F , 



n(kw) = -i[4(l- 7k )(J|G 1 |+<F 1 )]- 1 



A 2 = 



x / dte iut (A 2u Al), 
Jo 

i 2 si - w£s k . 



As follows from the above equation, to calculate w k 
and A2 we have to select terms of i 2 s k which are propor- 
tional to s k . It can be done only approximately because 
the quantit y J -k) cannot be calculated exactly. Fol- 
lowing Refs. ll4llfil we used the decoupling in i 2 s k for such 
selection and found 

col - 16aJ 2 (jGxl + (1 - 7 k)(A + 1 + 7k), (13) 

where A is the parameter of the gap in the spin exci- 
tation spectrum at the wave vector Q of the Brillouin 
zone. In an infinite 2D lattice this gap is opened for any 
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nonzero temperatureiS and at T = for x > 0.021^ The 
gap size is directly connected with the spin correlation 
length of the short-range antiferromagnetic order. Hence 
a finite gap for T > is in agreement with the Mermin- 
Wagner theorem*^ The gap parameter can be expressed 
through the model parameters and correlations of hole 
and spin operators^iS However, due to strong depen- 
dencies of the considered quantities on this parameter we 
found it more accurate to determine this parameter from 
the constraint of zero site magnetization (s£) = which 
is fulfilled in the paramagnetic state. This constraint can 
be written in the form 



V 1 - x ) = Jj E / ^coth(^) B(kw). 



(14) 



In Eq. i|13|) . the parameter a is introduced to improve 
somewhat the results obtained with the decoupling and 
to take into account vertex corrections. In earlier worksi& 
where the analogous correction were used for the Heisen- 
berg model this parameter was determined from the con- 
straint Q14p. Due to comparatively weak dependencies 
of the considered quantities on this parameter we found 
it more appropriate to set a — 1.802 — 0.802 tanh(lCte) 
and to use the constraint for the calculation of A, as 
mentioned above. The expression given for a takes into 
account its value obtained in Ref. [Tj for finite damping 
of spin excitations and the weakening of the vertex cor- 
rections with doping. 

When selecting terms of i 2 s£ which have to be included 
into A2 we omitted terms proportional to i 2 , being moti- 
vated by our earlier result and by the results of the spin- 
wave approximation 18 which indicate that n(kw) has to 
be proportional to t 2 . An additional argument to omit 
these terms is that a part of them contains multipliers 
of the type ca^amo- the mean values of which are 
zero. Other terms of this type and a part of terms pro- 
portional to t J contain the hole operators with opposite 
spins, aJno-Om'^CT, which also give zero on averaging and 
therefore were omitted. Terms which are proportional to 
J 2 and describe multiple spin-excitation scattering pro- 
cesses were not included into A2 either - in this article 
only the decay of the spin excitation into the fermion pair 
is considered. This process is described by the following 
terms: 



UJ 



A 2 



ffkkik 2 — 7k 



gkk 1 k 2 akio- a k+k 1 -k 2 , 

kik 2 <T 
1 



x(7k 2 -k! -7k x -7k+k x -k 2 +7k+kJ, 

where the line over the operators indicates that in cal- 
culating thermodynamic averages with A2 by factoriza- 
tion, terms containing couplings of hole operators from 
the same A2 have to be omitted, since such processes 
have already been included into ui^s^.. Substituting A2 
into the above definition of n(kw), neglecting the differ- 
ence between An and A 2 (i) and using the decoupling we 



get 



Imn(kcj) 



8^ 2 J 2 



1 - exp(w/T) 



iV 2 (l- 7k )(J|C 1 |+tF 1 ) 



E 

kik 2 



2 

.9kkik 2 



X [1 - + UJ\ — UJ2)B(k 2 Ld2) 

A(k 1 uj 1 )A(k + ki - k 2 , u> + uji - oj 2 ) 

L(k\UJ\a)L()i + ki - k 2 , uj + ui\ - uj 2 , a) 

(15) 



Equation (|13fl is supposed to give a good approximation 
for the real part of the frequency of spin excitations and 
therefore only the imaginary part of n(kcj) will be con- 
sidered below. Notice that ImII(k(jj) is negative, finite 
for uj = and even with respect to the change of the sign 

of UJ. 

As seen from Eq. ()15f) . Imil(kw) is finite for k — > 0, 
whereas vanishes in this limit. Therefore the spin 
Green's function 112fl has a purely imaginary, diffusive 
pole near the T point, in compliance with the result of 
the hydrodynamic theory. 19 In the general case proper- 
ties of spin excitations near the M point differ essen- 
tially from those near T. In the calculations of Ref. ^] 
for the former excitations the real parts of frequencies 
were larger than their imaginary parts due to the spin 
gap. However, it is worth noting that in this comparison 
only the decay into two fermions was considered as the 
source of damping. Another source of damping multi- 
ple spin-excitation scattering - was neglected. However, 
even in the case of overdamped excitations with k w Q 
their frequencies will have real components due to the 
spin gap. 

To simplify further calculations we take into account 
that in the considered underdoped case the spin spectral 
function B(ku>) is strongly peaked near Q for uj w wq. 
Allowing for the small value of u>q, A(k + Q, uj) w A(kuj) 
and Eq. ((Til) we get 



Imn(kw) 



97rf 2 J 2 (l - x) 



2N(l-^)(J\C 1 \ + tF 1 ) 
x E( 7k + k ' "7k') 2 



,n F (uj + uj') — tif(w') 
dw 

x A(k'w')4(k + k> + «') 
-L(k'u)'a)L(k + k',uj + u)',a)\. (16) 



Now the damping has taken the familiar form given by 
the fermion bubble. 
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IV. MAGNETIC SUSCEPTIBILITY 



We have used hole self-energies (J7J) and correlations of 
hole and spin operators obtained in Ref. ^| for calculat- 
ing the hole Green's function G, Eq. This function 
and the spin gap parameters A obtained in Ref. ^| have 
then been applied for the calculation of the spin Green's 
function determined by Eqs. lfl"2")). fl^l and |TH)l. This 
latter function is connected with the magnetic suscepti- 
bility by the relation 

X z (ktu) = -Af/ B D(kLo), 

where llb is the Bohr magneton. The self-energies of 
Ref.hJ | were calculated for a 20x20 lattice with the pa- 
rameters t — 0.5 eV, J = 0.1 eV which correspond 
to hole-doped cuprates 13 *SS and for the ranges of hole 
concentrations and temperatures < x < 0.16 and 
O.Oli w 58 K < T < 0.2i w 1200 K. For several hole con- 
centrations we have checked now that the self-energies 
calculated for T = O.Oli remain practically unchanged as 
the temperature decreases to T = 0.003i « 17 K. There- 
fore we can use these self-energies also for T < O.Oli. 
For temperatures close to zero the superconducting gap 
A s was set to O.OAt — 20 meV, the value extracted from 
the tunnelling experiments^ As follows from the experi- 
ments, this value remains practically unchanged with the 
doping variation from heavily underdoped to optimally 
doped YBa2Cu3C>7_y. 

Results of such calculations for the imaginary part of 
the magnetic susceptibility at the antiferromagnetic wave 
vector Imx(Q) are shown in Figs.^ an d|21 In these fig- 
ures experimental data™ on the magnetic susceptibility of 
underdoped YBa 2 Cu307_ y are also depicted. The oxy- 
gen deficiencies y — 0.5 and 0.17 in this crystal corre- 
spond to the hole concentrations x « 0.075 and 0.14^2, 
YBa2Cu307_j, is a bilayer crystal and the symmetry al- 
lows one to divide the susceptibility into odd and even 
parts. For the antiferromagnetic intrabilayer coupling 
the odd part can be compared with our calculations car- 
ried out for a single layer. 

The value of damping |Imn(Qw)| depends on widths 
of peaks in the hole spectral functions near the Fermi 
surface. These widths are determined by an artificial 
broadening which was introduced in Ref. to stabi- 
lize the iteration procedure. From the comparison with 
photoemission spectra™ of YBa2Cu30 , 7_ !/ it is seen that 
the peaks in Ref. ^| are more intensive and narrower 
than in experiment which leads to a larger value and 
stronger frequency dependence of the calculated damp- 
ing. To weaken this difference and to obtain a better fit 
of the shapes of the calculated susceptibility to the exper- 
imental data we have decreased |Imn(Qcj)| by a factor 
/ and added a constant damping r\ to it. This allows 
us to weaken somewhat the frequency dependence of the 
total damping T(Qo;) = |Imn(Qw)|//+?7. As will be dis- 
cussed in greater details later, the low-frequency shoulder 
in Imx(Qcij) is connected with this dependence. Thus, 
the fitting parameters / and r\ allow us to change the 



FIG. 1: The imaginary part of the spin susceptibility at the 
antiferromagnetic wave vector in the superconducting state. 
Curves show the results of our calculations in a 20x20 lat- 
tice for t = 0.5 eV, J = 0.1 eV, T = 17 K, x = 0.06 (a) 
and x — 0.12 (b). Filled squares are the odd susceptibility 
measured^ in YBa2Cu306.s (a, T c = 45 K, x w 0.075) and in 
YBa 2 Cu 3 06.83 (b, T c = 85 K, x w 0.14) at T = 5 K. Here 
and in Fig. |5|tick labels on the vertical axes correspond to the 
curves. In both figures experimental values are approximately 
1.5 times smaller than the calculated ones. 



FIG. 2: The imaginary part of the spin susceptibility in the 
normal state. Curves show the results of our calculations for 
T — 116 K, all other parameters are the same as for the 
respective panels in Fig. Q Filled squares are the odd suscep- 
tibility measured 2 in YBa2Cu30s.5 (a) and in YBa2Cu306.83 
(b) at T = 100 K. 
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FIG. 3: The frequency dependence of the total damping T(Q) 
used in the calculation of the four curves in Figs. □ and |2 
The parameter / is equal to 2.7 for the dashed curve and 
2 for the other curves. The parameter r\ is equal to 0.027t, 
0.04t, 0.012t, and 0.029* for the solid, dashed, short-dashed 
and dash-dotted curves, respectively. 



FIG. 4: The Fermi surface of the t-J model (lines) and the 
momenta which give the main contribution to the maxima of 
r(Q) in the used 20x20 lattice (circles). The antiferromag- 
netic wave vector Q connecting momenta of the fermion pair 
in the spin polarization bubble is shown by the arrow. The 
point M corresponds to k = (ir,ir). 



relative intensity of this shoulder. The damping 77 can be 
connected with the processes of multiple spin-excitation 
scattering or scattering at impurities. The frequency de- 
pendencies of the total damping used in the calculation 
of the curves in Figs. ^ and [21 are shown in Fig. [21 No- 
tice that the fitting parameters / and r\ with the values 
given in the caption to this figure influence only weakly 
the position of the maximum in susceptibility which is 
determined by the value of ujq. 

As seen from Figs.^ an d[21 the position of this maxi- 
mum, the resonance peak, and its evolution with doping 
and temperature described by the t-J model are in good 
agreement with those observed in YBa2Cu307_j,. In the 
model the maximum is connected with the excitation of 
localized Cu spins at the antiferromagnetic wave vector 
Q. Its frequency wq determines the size of the spin gap. 
In the underdoped case it determines also the frequency 
of the resonance peak. As shown in Ref. 0, wq grows 
with doping and this leads to the growth of the frequency 
of the resonance peak from approximately 18 meV at 
x = 0.06 to 38 meV at x = 0.12 in Figs. □ and It 
was also showni4 that Imx(kw) is strongly peaked at Q 
and that the value of Imx(Qw) decreases with doping 
which is in agreement with experimental observations. 1 '- 2 
In absolute units our calculated values of Imx(Qw) are 
approximately 1.5 times larger than its experimental val- 
ues. 

We notice that the shape of the calculated frequency 
dependence of the susceptibility is close to that observed 
experimentally. Of special interest is the low-frequency 
shoulder in this dependence. This shoulder is more pro- 
nounced for lower hole concentrations and temperatures. 
As mentioned above, it originates from the strong fre- 
quency dependence of the damping T(Q) shown in Fig. [3 



The pronounced maxima of the curves in this figure 
are connected with intensive peaks in the hole spectral 
function for momenta near the Fermi surface. These 
peaks correspond to the so-called spin-polaron bandit 
For moderate doping the Fermi surface of the t-J model 
consists of two rhombuses with rounded corners^ These 
rhombuses are centered at the T and M points and are 
approximately nested by the momentum Q. This nesting 
is also very essential for the appearance of the maximum 
in r(Qw). In Fig.0]the Fermi surface is shown and mo- 
menta of the hole spectral functions which give the main 
contribution to the maxima of T(Q) in the used 20x20 
lattice are indicated. For these momenta the intensive 
spin-polaron maxima in the spectral functions AQc'u)') 
and A(Q + k', u> + a/) [see Eq. (I1G[) ] overlap and fall into 
the frequency window determined by the difference of the 
occupation numbers. 

The Fermi surface in YBa2Cu307_ a differs from that 
shown in Fig. ^J 23 - 24 However, it is known from the pho- 
toemission experiments that at least in the superconduct- 
ing state the hole spectral function has pronounced peaks 
for momenta near the Fermi surface. In the two-layer 
YBa2Cu307_ a the main contribution to the damping of 
the spin excitations is given by the decay into the fermion 
pair in which one of the fermions belongs to the bond- 
ing band and the other to the antibonding band and the 
respective parts of the Fermi surface are nested by the 
momentum (w, it, it) These conditions are similar to 
those observed in the t-J model and therefore the low- 
frequency shoulder in the susceptibility in YBa2Cu307_ y 
can be also related to the strong frequency dependence of 
the damping of the spin excitations which arises due to 
pronounced peaks in the hole spectral function and the 
nesting. 
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FIG. 5: Constant energy (k, k) scans at the resonance energy 
ujq. Solid and dashed curves show the results of our calcula- 
tions for x = 0.12, ujq — 38 meV in the superconducting state 
at T — 17 K and in the normal state at T = 116 K, respec- 
tively. To simulate a finite instrumental momentum resolu- 
tion the curves were calculated by the convolution of Imx(ko;) 
with the Gaussian with the full width at half maximum equal 
to 0.2-k in the momentum space. Filled and open squares 
are experimental data— in YBa2Cu3 06.83 for wq = 35 meV, 
T = 4 K and 109 K, respectively. 

It is worth noting that for all four curves in Figs. ^ 
and |21 the value of T(Q, cjq)/2 is smaller than ujq. Thus, 
in contrast to a vicinity of the T point near the M point 
the spin excitations are not overdamped in underdoped 
YBa 2 Cu 3 7 - y . 

Now let us consider the momentum dependence of 
the resonance mode. In Fig. [3] the constant energy 
scans obtained in our calculations are compared with 
experiment in YBa 2 Cu306.83- The scans were per- 
formed along the diagonal of the Brillouin zone at the res- 
onance energy in the superconducting and normal states. 
To simulate a finite instrumental momentum resolution, 
which is comparable to the width of the peak in Imx(kcj) 
our curves were calculated by the convolution of this 
quantity with the Gaussian with the full width at half 
maximum equal to 0.2ir in the momentum space. This 
corresponds to 0.1 in reciprocal lattice units which is the 
usual resolution in experiments of this type. As can be 
seen from Fig. [H] for both temperatures the calculated 
momentum dependencies are in good agreement with ex- 
periment. 

In Fig. the dispersion of the maximum of our calcu- 
lated susceptibility is compared with experimental dataS 
in YBa2Cu3C>6.5. This dispersion corresponds approx- 
imately to Wk in Eq. i|13|) ■ Fot small q = k Q this 
momentum dependence can be written as 

Wk » ^ + C 2(k-Q)2. (i 7 ) 

This function fitted to our calculated data with the pa- 
rameters ujq = 18.4 meV and c/a — y/8a\Ci\J = 



FIG. 6: The dispersion of the maximum in the frequency 
dependence of Imx(qo;), q = k — Q. Filled squares are our 
results for x — 0.06 and T = 17 K. The fit for these data 
with Eq. 1171 is shown by the curve. Open squares are ex- 
perimental results- in YBa2Cu30e.5 at T = 5 K for odd spin 
excitations. 

0.134 meV is also shown in Fig. |SJ Here a is the dis- 
tance between Cu sites in a Cu-0 plane. As seen from 
this figure, our calculated dispersion is close to the exper- 
imental one for similar parameters. For u> « ujq Imx(kaj) 
is peaked at k = Q. For to > ujq the susceptibility has 
maxima on the ring with the radius approximately deter- 
mined by the equation u — [wq + c 2 (k — Q) 2 ] 1 / 2 . In con- 
stant energy scans along some direction this property of 
the susceptibility manifests itself as two peaks in incom- 
mensurate positions equally spaced from the M point. 26 
For lj < loq for the considered parameters Imx(ku) is 
peaked at k = Q. 



V. CONCLUDING REMARKS 

We have considered the magnetic susceptibility for the 
underdoped case when the resonance peak is observed 
both in the normal and in the superconducting states. 
As mentioned, the frequency of the peak is determined 
by the frequency of the spin excitation ujq which sets the 
size of the spin gap. This frequency grows with the hole 
concentration^ in agreement with experimental obser- 
vations in underdoped crystalsi 2 * 6 - 

For the normal-state t-J model in the overdoped region 
it was shownS that the part of the magnon branch, which 
persisted at lower doping at the periphery of the Brillouin 
zone, is suddenly destroyed for x w 0.17 at T = 0. This 
transition is accompanied by the radical change of the 
hole spectrum: dispersion and distribution of the spec- 
tral weight become close to the case of weakly correlated 
fermions. This result corresponds to the sudden disap- 
pearance of the resonance peak in the normal-state over- 
doped cuprates. 2,6 One of the reasons for the transition 
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in the t-J model is the damping of the spin excitations 
which grows with doping. A considerable decrease of the 
damping in the superconducting state can restore the 
spin excitations near the M point in the frequency range 
lu < 2A S . Such mechanism was considered in Ref. [ll| 
where the magnetic susceptibility similar to that given 
by Eqs. I|12|) and (|13|l was postulated and the damping 
described by the fermion bubble of the type of Eq. i|16|) 
was used. Above the mentioned transition at x ~ 0.17 
the hole spectrum of the t-J model becomes similar to 
that used in Ref. O an d the analogous outcome can be 
expected here. 

In contrast to the underdoped region, in the overdoped 
case the frequency of the resonance peak decreases with 
doping which can be related to a finite damping of the 
spin excitations and to the decrease of the superconduct- 
ing gap with doping in this range of concentrations 

In contrast to YBa2Cu307_ y where ojq < 2A S , in 
La 2 _ 2; Sr :I .Cu04 the value of 2A S « 9 meV is substantially 
smaller than wq which is supposed to be approximately 
the same as in the former crystal. This difference may 
be the reason for the absence of the resonance peak in 
overdoped L^-xSr^CuC^iii Changes in the susceptibil- 
ity observed^ in Lai.86Sr .i4CuO4 at the superconduct- 
ing transition consist of some suppression of Im% below 
the superconducting gap and an increase above it. The 
suppression can be connected with the decrease of the 
damping of the spin excitation accompanying the open- 



ing of the gap, while the increase of the signal above the 
gap is apparently a combined effect of the transfer of 
the carrier spectral weight above the gap and the nesting 
supposed— for the Fermi surface of this crystal. 

In conclusion, we compared the magnetic susceptibility 
calculated in the t-J model with the experimental data 
in the underdoped YBa2Cu307_j,. It was demonstrated 
that the calculations reproduce correctly the frequency 
and momentum dependencies of the experimental sus- 
ceptibility and its variation with doping and temperature 
in the normal and superconducting states. This allowed 
us to interpret the maximum in the frequency depen- 
dence - the resonance peak - as a manifestation of the 
excitation branch of localized Cu spins and to relate the 
frequency of the maximum to the size of the spin gap. 
The low-frequency shoulder well resolved in the suscep- 
tibility of superconducting crystals was connected with a 
pronounced maximum in the damping of the spin excita- 
tions. This maximum is caused by intense quasiparticle 
peaks in the hole spectral function for momenta near the 
Fermi surface and by the nesting. 
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